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Exact calculated time evolutions in the framework of a many-electron model of itinerant mag¬ 
netism provide new insights into the laser-induced ultrafast demagnetization observed in ferromag¬ 
netic (FM) transition metal thin films. The interplay between local spin-orbit interactions and 
interatomic hopping is shown to be at the origin of the observed post-excitation breakdown of 
FM correlations between highly stable local magnetic moments. The mechanism behind spin- and 
angular-momentum transfer is revealed from a microscopic perspective by rigorously complying with 
all fundamental conservation laws. An energy-resolved analysis of the time evolution shows that the 
efficiency of the demagnetization process reaches almost 100% in the excited states. 
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Pump-and-probe femtosecond-laser experiments on 
thin magnetic transition metal (TM) films have shown, 
almost 2 decades ago, that the ferromagnetic (FM) or¬ 
der breaks down on a time scale of a few hundred fem¬ 
toseconds after the pulse absorption.® This remarkable 
finding opened the way to an ever since growing re¬ 
search field, which has a wide fundamental and practi¬ 
cal importance.lSHIll The phenomenon as such can be nei¬ 
ther an immediate consequence of the excitation, since 
optical transitions conserve spin, nor the result of ther¬ 
mally activated stochastic processes, which would involve 
a much longer time scale. Instead, ultrafast demagnetiza¬ 
tion (UFD) reflects the intrinsic many-body dynamics of 
correlated excited electrons in FM metals. Understand¬ 
ing its origin and controlling its properties is therefore of 
crucial importance.li^ 

From a fundamental perspective, it has been clear from 
the start that spin-orbit (SO) interactions must play a 
central role in the dynamics, since only the relativistic 
corr ections break the conservation of the total electronic 
spinLilliiTs] jndeed, the spin-orbit coupling (SOC) allows 
exchanges between the dominant local Zd spin moments 
Si and the local orbital moments U at every TM atom 
i. Taking into account that the total angular momentum 
ji = h + Si is conserved in this process, the focus of atten¬ 
tion quickly moved towards quantifying the time depen¬ 
dences of the spin and orbital moments, and to correlate 
them to the observed demagnetization. This has been 
experimentally achieved by performing time-resolved X- 
ray magneti c circular dichroism (XMCD) measurements 
on These works showed, in contrast to early ex¬ 

pectations, that no enhancement of liz accompanies the 
decrease of Siz, but rather that both Siz and Uz decrease 
as a function of time with a time constant of r ~ 120 fs. 
Since li is not a reservoir for angular momentum, the 
authors concluded that a femtosecond spin-lattice relax¬ 
ation, i.e., a substantial femtosecond spin angular mo¬ 
mentum transfer to the lattice, takes place.® 

Although the theoretical research in this field has been 
most intense, understanding the microscopic mechanisms 


of UFD and angular-momentum relaxation still remains 
an open problem.l^ Over the past years, two different 
theoretical approaches have attracted particular atten¬ 
tion. One of them is electron-phonon spin-flip scattering, 
in which the lattice is assumed to be a perfect sink for 
angular momentum.^iSHlIl The other one is spin-polarized 
electron diffusion, which does not invoke any angular mo¬ 
mentum dissipation channel, but rather a spin-dependent 
superdiffusive electron transport from the laser-excited 
film to the substrate.^^^^^ In this context it is quite re¬ 
markable that none of these theories happens to bear 
a clear relation with the fundamentals of itinerant mag¬ 
netism, which are tightly anchored to strong electron cor¬ 
relations and to the resulting high stability of the local Zd 
magnetic moments.!^ Only recently the potential impor¬ 
tance of l ocaliz ed spins and their fluctuations has been 
suggest ed.l22lSni it is therefore most challenging to estab¬ 
lish the links between equilibrium and non-equilibrium 
theories of itinerant magnetism. The purpose of this 
Letter is to develop a many-body theory of UFD, which 
takes into account the electronic correlations responsible 
for local moment formation and magnetic order, and to 
analyze its physical consequences rigorously by perform¬ 
ing exact time propagations. 

Several distinct features are expected to be central to 
the physics of the laser-excited electrons in FM met¬ 
als: (i) the single-particle hybridizations responsible for 
electron delocalization, bonding and metallic behavior, 
(ii) the dominant Coulomb interactions among the Zd 
electrons, which introduce correlations, stabilize the local 
magnetic moments and, together with the single-particle 
contributions, define the magnetic order, (iii) the SO in¬ 
teractions, which couple the spin and orbital degrees of 
freedom, and (iv) the interaction with the external laser 
field, which triggers the initial electronic excitation. We 
therefore propose the pd-band model given by 


H = H^ + Hc + Hso + HE{t) , 


( 1 ) 
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where 

Hq = Safliaa + tif (^) 

iQ;(T Zjf Ct0(T 

describes the band structure of the relevant 3d and Ap 
valence electrons that are responsible for the magnetic 
properties in TMs and for the relevant optical absorp¬ 
tion. The operator {ciaa) creates (annihilates) a 

spin-CT electron at atom i in the orbital a, which has 
well-defined radial and orbital quantum numbers nlm. 
fiiaa = ^iaa-^iacr counts the Corresponding occupations. 
The energy of the orbital a is denoted by £„ and the in¬ 
teratomic hopping integrals by The second term in 
Eq. 0 designates the dominant intra-atomic Coulomb 
interaction among the 3d electrons 

He = -1) - J^sf ■ sf , (3) 

i i 

where U stands for the average d-electro n dire ct Coulomb 
integral and J for the exchange integral.l^^^ The opera¬ 
tors nf and refer, respectively, to the d-electron num¬ 
ber and total spin at atom i. The spin-orbit interactions 
are given by the third term 

Hso = C , (4) 

i cx^G3d gg' 

where (Z • stands for the intra-atomic matrix ele¬ 

ments of Z • s and ^ is the SOC constant. For simplicity, 
Ap electrons are here ignored. Finally, the last term 

HE{t)=eE{t)-'^{a\?\l3)c\^^Cipa ( 5 ) 

ia^G 

describes the interaction with the external laser field E{t) 
in the intra-atomic dipole approximation (e > 0 is the 
electron charge). The usual atomic selection rules for the 
position operator f imply that only dp and pd transitions 
enter the sum. 

At this point it is useful to recall the fundamental 
conservation laws underlying the model, which are the 
same as in the first-principles Hamiltonian. The non- 
relativistic terms Hq, Hq and He conserve the total spin 
^ = Y.ta since [Hq, 5 ] = [He, 5 ] = [He, S] = 0 . The 
spin conservation is broken only by Hso oc + 

h-Si+)/2 + hzSiz], which involves intra-atomic angular- 
momentum transfer between Si and U. Still, the total 
angular momentum ji = h + Si is locally conserved in 
the SO transitions, since [Hso, Si -I- Z^] = 0 for any atom 
i. Intra-atomic Coulomb interactions are also invariant 
upon rotations and thus preserve s), U and ji. How¬ 
ever, the interatomic hybridizations, though total-spin 
conserving, do not preserve the orbital angular momen¬ 
tum L = since the hopping integrals connect 

orbitals having different m at different atoms {[Hq, T] 0 


and [Ho,li] 7^ 0).l^We shall see that these simple argu¬ 
ments allow us to understand a number of qualitative 
aspects of the magnetization dynamics, its dependence 
on the relevant physical parameters, and the main mech¬ 
anism behind it. 

Despite the simplicity and transparency of the pro¬ 
posed peZ-band model, an analytical or straightforward 
numerical solution of its dynamics is out of reach at 
present. As in the equilibrium case, the main difficul¬ 
ties stem from the the Coulomb interaction Hq and the 
resulting many-body behavior. One could in principle re¬ 
sort to time-dependent mean-field approximations. How¬ 
ever, these are known to introduce artificial symmetry 
breakings, which spoil the fundamental spin rotational 
invariance, thus casting potentially serious doubts on the 
resulting magnetization dynamics. In order to derive rig¬ 
orous conclusions, we have therefore decided to consider 
a simplified small-cluster version of the model and to ob¬ 
tain an exact numerical solution of the ground state, ex¬ 
citations and time propagation. Similar approaches have 
been most successful in the context of equilibrium proper¬ 
ties of narrow-band systems.l^^^^^In addition, our results 
show that the physics of the magnetization dynamics can 
be explained by short-range effects, which justifies the 
small-cluster approach a posteriori. 

The parameters used for the calculations correspond 
approximately to Ni and have been specified as follows. 
The nearest-neighbor (NN) Slater-Koster hopping inte- 
grals are taken from band structure calculations .1^ The 
dp promotion energy Acp^ = 1 eV yields a dominant 
3cZ band occupation in the ground state. The direct 
Coulomb integral U = 4.5 eV and the exchange integral 
J = 0.8 eV lead to FM order with nearly saturated 3cZ 
moments at low temperatures.!^ Finally, the SOC con¬ 
stant ^ is varied in the range |^| < 100 meV, which in¬ 
cludes the values found in 3cZ TMs.!^ For the numeri¬ 
cal calculations we reduce the degeneracy of the bands, 
by considering only three 3cZ orbitals per atom (having 
m = — 1, 0 and 1) and one Ap orbital per atom (m = 0). 
The fee (111) monolayer is modeled by an equilateral 
triangle having A^e = 4 valence electrons. Since the aver¬ 
age occupation of the d orbitals is below half-band filling 
{{nfa) — He/iNa = 0.44) we set ^ < 0 in order that Si 
and li align parallel to each other.l^ 

A first test on the validity of the model and parame¬ 
ter choice is provided by the ground-state results, which 
match qualitatively the magnetic properties of Ni. From 
exact Lanezos diagonalizations we obtain that the ground 
state is FM with an off-plane easy magnetization axis and 
a magnetic anisotropy energy AE = 2.8 meV per atom. 
The local spin momenta 2siz = 1.32H are almost satu¬ 
rated, whereas the local orbital momenta kz = 0.09Zi are 
quenched to a large extent. These values should be com¬ 
pared, for example, with 2siz = 0.62h and Z^j. = 0.07Zi 

obtained in experiment.^^S! 

Starting from the FM ground state we determine the 
exact time evolution numerically by performing short- 
time iterative Lanezos propagations.^!! The actual dy- 
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FIG. 1. (Color online) Time dependence of the spin mag¬ 
netization 2Sz following a 5 fs laser-pulse excitation for an 
equilateral triangle having A^e = 4 valence electrons and dif¬ 
ferent SOC strengths The inset shows the demagnetization 
time Tdm as a function of 

namics is triggered by the femtosecond laser pulse E{t) = 
e ■ ilo cos(a;t) exp(—whose polarization vector e 
is along a NN bond within the (111) plane. The laser 
wave length A = 800 nm corresponds to the p hoton en- 
ergy hui = 1.55 eV used in several experiments.l^EISIMIini 
The pulse is centered at t = 0 and has a duration of 
Tp = 5 fs. Its amplitude Eq is such that the electronic 
system absorbs 0.2 photons per atom on average, which 
corresponds to a fluence F ~ 40 mJ/cm^ 

In Fig. results are given for the magnetization dy¬ 
namics. For realistic values of |^| = 50-100 meV the 
spin angular momentum per atom relaxes irreversibly 
from 2Sz/Na = l.32h to around 0.641i within the first 
hundreds of femtoseconds following the laser excitation. 
This demonstrates the ultrafast demagnetization effect in 
agreement with experiment.^ Fig. also reveals the cen¬ 
tral role played by the SOC. For ^ = 0 the total spin is un¬ 
affected by the optical excitation, as expected. Moreover, 
as the SOC is turned on one observes that the degree 
of demagnetization ASz = Sz{0) — ^^(oo) increases and 
that the demagnetization time Tdm decreases. Finally, for 
larger |^| > 50 meV the SOC strength affects only the 
time scale. In order to quantify Tdm, we have fitted the 
exact time dependence Sz(t) with an exponential func¬ 
tion of the form (t) = exp(—t/rdm) + <5'z(0) — 
which is shown by the dashed curves in Fig. The ob¬ 
tained Tdm, shown in the inset, behaves approximately 
as Tdm ~ 6/i/|^|. This is consistent with the time-energy 
uncertainty relation and confirms that the SOC, which 
represents the smallest energy scale, controls the time 
scale of the relaxation process. 

The time dependences of the local d-electron spin mo¬ 
ments = yj ((s/)^) and of the NN spin correlation 
functions shown in Fig.^for ^ = —50 meV 

provide a much clearer picture of how xhe UFD actually 
occurs. One observes that fj-i changes very slightly from 
IMh to 0.98fi, while the laser pulse is on (r^ = 5 fs). This 
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FIG. 2. (Golor online) Time dependence of the average d- 
electron local spin moment ((s)”^)^)^^^, NN spin-spin correla¬ 
tion function (s/ • s/), and off-plane spin and orbital angular 
momentum Sz and Lz. 


is the result of a small laser-induced dp charge transfer, 
which involves majority d electrons, and which causes 
the number of d-electrons per atom to decrease from 
Ud ~ 1.25 to Ud — 1.17. Otherwise, /if is essentially 
time independent. The remarkable stability of the lo¬ 
cal d moments shows, as in thermal equilibrium above 
the Curie temperature Tc, that the laser-induced ultra¬ 
fast breakdown of FM order is not the consequence of a 
significant loss of local spin polarization.^^ Instead, the 
UFD follows to a large extent from the quantum fluctu¬ 
ations of the orientations of the local magnetic moments 
sf, whose magnitude remains almost constant through¬ 
out the process. Indeed, as shown in Fig. the average 
NN spin-correlation function decreases very rapidly 
as a function of time, approaching its long-time limit al¬ 
ready 50-80 fs after the laser absorption, and remaining 
approximately constant in the following. The breakdown 
of the NN correlations is actually much more rapid than 
the decrease of the total spin polarization Sz (see Fig.[^. 
The same trend holds for other values of for example 
^ = —100 meV, where the demagnetization time is much 
shorter. This shows that the quantum fluctuations of sf 
dominate only the first stages of the spin dynamics. A 
slower tilting of S off the z axis follows, during which the 
transversal components Sx and Sy always remain zero. 

In order to analyze the mechanisms of angular- 
momentum transfer behind the UFD process, we turn 
our attention to the dynamics of the orbital moment Lz 
and contrast it with the dynamics of Sz (see Fig. 1^. In 
the ground state, before the laser excitation, Lz/Na — 
0.09h is quenched to a large extent in comparison with 
the Hund-rule atomic value Lf- = IH, as expected in 
TMs. After the excitation Lz/Na shows rapid oscilla¬ 
tions between 0.02h and 0.12h, always remaining parallel 
to Sz- This oscillatory behavior is numerically stable 
and perfectly reproducible. The total angular momen¬ 
tum Jz = Lz + Sz is obviously not a constant of motion. 
This is in agreement with time-resolved XMCD experi¬ 
ments showing that Lz is not a reservoir for the decreas- 
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FIG. 3. (Color online) Time dependence of the spin mag¬ 
netization 25^“^ in the excited-state manifolds corresponding 
to a = 0, 1 and 2 photon absorptions {Tilj = 1.55 eV). The 
average Sz = Tla is given by the dashed curve. The 

inset shows the spectral distribution Dis,{e) of 4^(1) after the 
laser-pulse passage {t > Srp = 15 fs). The weights Wa of the 
different spectral parts of 'F are indicated. 

ing Sz 2^2 remains quenched for all 

times is the result of the interatomic hybridizations re¬ 
sponsible for the electronic motion in the lattice and for 
the band formation. Formally, one could say that the op¬ 
erator Ho preserves S but not L (i.e., = 0 while 

[Hq, L] yf 0) since the potential generated by the ions is 
not rotationally invariant But physically it is impor¬ 
tant to realize that the characteristic time Tg required 
to quench the orbital moment of an electron in a lattice 
is extremely short, of the order of Tq ~ fi/Wd — 0.1 fs 
where Wd — 6 eV stands for the d-band width.SS The 
combination of a local J-conserving transfer of angular 
momentum from S to L, due to SO interactions, and 
a very fast dynamical quenching of L, due to electronic 
hopping, explains the spin-to-lattice angular-momentum 
transfer observed in experiment. In the framework of 
the present model, the dynamical quenching of L would 
manifest itself only as a global rotation of the rigid lat¬ 
tice, since electron-phonon coupling (EPC) has been ne¬ 
glected. Including EPC in the model would open an ad¬ 
ditional channel for L-quenching, through which angular 
momentum would be transferred to the lattice vibrations. 

In order to verify the validity of our conclusions we 
have repeated the calculations of the dynamics by re¬ 
ducing artificially the hopping integrals and thus ap¬ 
proaching the atomic zero-band-width limit. As shown 
in the Supplemental Material, one observes that as soon 
as becomes comparable or smaller than the SOC con¬ 
stant ^ the rapid quenching of is replaced by oscilla¬ 
tions of both Lz and Sz, with a period of about 70 fs, 
during which Jz is approximately conserved.!^ 

A complementary perspective to the UFD process is 
obtained by performing a spectral analysis of the laser- 
excited many-body state 4>(t) and of the time depen¬ 


dence of Sz- For this purpose, we considered 41(t) imme¬ 
diately after the excitation (at t = Srp) and expanded it 
in the stationary states \tjjk) of the field-free Hamiltonian 
H = Hq + He + Hso having energy e^. The obtained 
spectral distribution D^{e) = ~ ^fc)K''/'fc I ^ )P 

is shown in the inset of Fig. Notice that D^,{e) is 
independent of t after the pulse passage, since so is H 
[see Eqs. 0”(§]- Three main groups or manifolds of 
nearby peaks can be recognized, which are located at the 
excitation energies As ~ ahuj and which correspond to 
the absorption of a = 0, 1 and 2 photons. The spin 
magnetization si°‘'^ originating from these manifolds is 
shown in Fig. as a function of t. One observes that 
immediately after excitation all a have the same satu¬ 
rated magnetization. This simply reflects the spin con¬ 
servation upon optical dipole transitions. The ground 
state, being a pure stationary state, preserves {Sz) and 
yields a time-independent si^'^. Relaxation can only stem 
from the excitations. Indeed, in the excited manifolds 
si°‘'^ decreases dramatically to nearly zero, particularly 
in the most relevant a = 1 manifold. This demonstrates 
the remarkably high efficiency of laser-induced UFD in 
the excited states. The residual average magnetization 
2Sz{t —)■ oo), which persists after the spin-lattice relax¬ 
ation process (see Fig. is essentially the consequence 
of the finite overlap between the many-body state 4'(t) 
after the laser absorption and the ground state. 

In sum, the laser-induced magnetization dynamics of 
ferromagnetic TMs has been studied in the framework of 
an electronic model. For the first time a solution of the 
time-dependent many-body problem has been achieved, 
which is complete from the perspectives of electron cor¬ 
relations, spin-orbit interactions and essential symme¬ 
tries. The results have demonstrated that the femtosec¬ 
ond demagnetization can be explained in terms of a three- 
step mechanism: (i) The laser pulse creates electron-hole 
pairs. This opens the way for (ii) the SOC yielding local 
angular-momentum transfer from Si to U with a charac¬ 
teristic time scale of fi/\^\ « 10 fs. However, angular mo¬ 
mentum is not accumulated in U, since (iii) L is quenched 
by the motion of electrons in the lattice. This takes place 
on a much shorter time scale of only h/Wd ^ 1 fs. Exten¬ 
sions of this work by improving the model Hamiltonian 
are certainly desirable and necessary for a more realis¬ 
tic description of specific magnetic materials. In par¬ 
ticular the coexistence of different relaxation channels 
(electronic relaxation, electron-phonon spin-flip scatter¬ 
ing, spin diffusion, etc.) deserves to be explored. In any 
case, the simplicity of the electronic processes identified 
in this work and the fundamental character of the pro¬ 
posed model suggest that the rigorously derived concepts 
should be universally applicable. 
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